Superfluid vortex front at T — )■ 0: Decoupling from the reference frame 
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Steady-state turbulent motion is created in superfluid ^He-B at low temperatures in the form of 
a turbulent vortex front, which moves axially along a rotating cylindrical container of '^He-B and 
replaces vortex-free flow with vortex lines at constant density. We present the first measurements on 
the thermal signal from dissipation as a function of time, recorded at 0.2 Tc during the front motion, 
which is monitored using NMR techniques. Both the measurements and the numerical calculations 
of the vortex dynamics show that at low temperatures the density of the propagating vortices falls 
well below the equilibrium value, i.e. the superfluid rotates at a smaller angular velocity than the 
container. This is the first evidence for the decoupling of the superfluid from the container reference 
frame in the zero-temperature limit. 

PACS numbers: 67.30. hb, 02.70.Pt, 47.15.ki, 67.30. he 



Introduction:-For eighty years the phenomenologi- 
cal two-fluid model of Tisza and Landau has been the 
starting point for discussions on dynamics in superfluids. 
Onsager and Feynman introduced the quantized vortex, 
which plays an all important role in the dynamics: Vor- 
tices mediate the interaction between the superfluid and 
normal components of the liquid via reactive and dissipa- 
tive friction forces [ij. Vortices also allow the superfluid 
component to participate in non-potential flow and thus 
in rotation. 

The problem of rotation is important in modern 
physics m. It goes back at least to Newton and his 
bucket experiment: In a container which rotates at con- 
stant angular velocity f2 with respect to the inertial frame 
a viscous liquid finally relaxes to the equilibrium state - 
the state of solid-body rotation with the container. This 
state is distinguished eg. by the parabolic shape of the 
meniscus, which allows to select the absolute reference 
frame in which the surface of the liquid is flat. The su- 
perfluid provides new insights to the problem of rota- 
tion: At r = (where the normal component is absent) 
a vortex-free superfluid remains at rest in the laboratory 
frame, its meniscus is flat and thus does not discriminate 
between different rotating frames. As superfluids have 
been found to simulate in some cases the quantum vac- 
uum the modern ether 0, , one can ask is rotation an 
absolute effect in the quantum vacuum 0] . This reminds 
us about the arguments of Mach, who claimed that rota- 
tion is a relative effect and can be distinguished only due 
to coupling to another reference frame (distant stars). 

Quantized vorticity changes the situation: by acquir- 
ing an array of rectilinear vortex lines the superfluid com- 
ponent is able to participate in the rotational flow and 
to mimic solid-body rotation. The friction force between 
vortices and the normal component provides the coupling 
between the container frame and that of the rotating su- 
perfluid, so that in complete equilibrium these frames co- 



incide: vortices arc at rest in the container frame. How- 
ever, the coupling between the frames decreases with de- 
creasing temperature, and is lost in the T — >■ limit 
when the normal component vanishes. Here we discuss 
an experimental example of the decoupling of the frames 
- the propagating vortex front in superfluid ■^He-B. We 
find that in the low-T limit, the superfluid component be- 
hind the front develops its own rotating reference frame, 
whose angular velocity is smaller than that of the con- 
tainer and is decreasing with decreasing coupling. 

Turbulent vortex-front motion is ubiquitous in the dy- 
namics of both viscous fluids Q and superfluids . In 
superfluid '*He turbulent fronts and plugs are known from 
pipe flow driven by thermal counter currents of the nor- 
mal and superfluid components. In a long cylinder of 
superfluid "^He-B the spin-up response of the superfluid 
component from rest to rotation often takes place as an 
axially propagating precessing vortex front |7[. The mo- 
tion of the front is launched after a sudden localized tur- 
bulent burst. This can be engineered to happen at con- 
stant angular velocity fl, if the superfluid is initially in 
the metastable vortex-free Landau state and the turbu- 
lent burst is triggered externally. Once formed, the front 
consists of a turbulent core, with an axial length compa- 
rable to the cylinder radius R, and a quasi-laminar tail in 
the form of a helically twisted bundle Q of vortex lines. 
The front is one of the few examples of steady-state tur- 
bulent motion in superfluids, which are presently avail- 
able for experimental investigations in the T — > limit. 
We have performed the flrst direct measurement of the 
heat released by such quantum-turbulent motion. 

The difference in the free energy of the vortex-free su- 
perfluid and that in equilibrium rotation is npsR^fl^ /4: 
per unit length of the cylinder of radius R. Here ps is 
the density of the superfluid component and a small cor- 
rection from the discrete nature of vortices is neglected. 
In our measurements the spin-up process happens at 
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FIG. 1: Measuring setup. Two independent cw NMR spec- 
trometers are used to monitor vortex front motion in the sam- 
ple container above the 0.75 mm orifice. The middle section 
with the quartz tuning fork oscillators is employed for ther- 
mometry and quasiparticle bolometry. The bottom section 
below the 0.3 mm orifice provides the thermal contact to 
the refrigerator. The superconducting solenoid in the center 
of the NMR sample produces the field for stabilizing the 
^He-A phase. This A-phase barrier layer divides the NMR 
sample in two "^He-B sections of equal length. 



f2 = const while the normal component is corotating with 
the container. In this case the work performed by the ro- 
tation drive can be found using the balance of angular 
momentum transfer. It can then be shown that the heat, 
expected to be dissipated in the spin-up process, is equal 
to the change of the energy of the superfluid. 

The simplest model of stationary-state front motion is 
that of a thin turbulent front to which all dissipation is 
concentrated and where the laminar relaxation to recti- 
linear vortex lines behind the front is neglected. If the 
front moves with the velocity V ^ the dissipation rate is 
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For stationary state motion Eq. ([T]) gives the maximum 
possible rate of heat release. At 0.25 and below the 
front velocity was found in Ref. to be given by F ss 
Q!cff rji?, where Ofoff ^ 0.1 is a constant, generated by the 
turbulence in the front. In this case the signal intensity 
\& Q with only weak temperature dependence. 

Experimental techniques:-In Fig.[T]the ■^He-B sam- 
ple in the top section of the long sample cylinder can 
be organized to rotate around its symmetry axis at con- 
stant angular velocity i7 in the vortex- free state Q . Here 
the superfluid fraction is at rest in the laboratory frame 
(vg = 0), while the normal excitations are in solid-body 
rotation (v,i = il x r). The distribution of vortices and 
superfluid counterflow as well as the motion of vortices 
is surveyed with two NMR detector coils • To trigger 
the front propagation, we use the Kelvin- Helmholtz shear 
flow instability of the AB interface [ll|. The instabil- 
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FIG. 2: Two fi-scaled bolometer signals. Time f = is as- 
signed to the moment when the trigger signal (cooling due to 
the latent heat in B to A transition) is observed. This trig- 
ger signal and the background heat leak has been removed by 
subtracting the equivalent signal measured in the equilibrium 
vortex state. The thermal time constant of the bolometer is 
about 25 s, which is much shorter than the time scale of the 
front propagation. 



ity occurs at a well-defined critical velocity ^ab, which 
depends on the magnetic stabilization field iJb and its 
gradient at the location of the AB interface. In our ther- 
mal measurements f^AB is traversed by sweeping _ffb at 
57 = const, for example, by increasing i/b from zero until 
the A phase is formed. The important consequence from 
the instability event is the escape of a number of vortex 
loops across the AB interface from the A to the B phase 
side. The loops interact and produce a large number of 
vortices in a turbulent burst close to the AB interface. 
In the setup of Fig. [T] the instability occurs simultane- 
ously at the two AB interfaces. Thus both upward and 
downward propagating fronts are set into motion. 

Thermal measurements:-The heat released in the 
spin-up process is recorded as a temperature rise across 
the lower orifice, which is the dominant thermal resis- 
tance in the ballistic regime. The value of this resistance 
is calibrated from the temperature increase as a function 
of the known heating power, using one quartz tuning fork 
oscillator as thermometer and the other as heater [T^ . 
The calibration also gives a residual heat leak to the '^He 
sample of 15 pW at 51 = 0, which increases by ^ 4pW 
for a Irad/s increment in fi. 

Examples of the thermal measurements of the front 
propagation are shown in Fig. [2] The total energy re- 
lease, integrated above background, matches the expec- 
tations: 0.55 and 0.75 nJ at 1 and 1.2rad/s, respectively, 
versus the expected values of 0.73 and 1.01 nJ. The dif- 
ference of 25% between the measured and expected heat 
realease can be ascribed to uncertainties in the bolometer 
calibration. 

The arrival of the front to the end of the sample, as 
determined from the NMR measurements, corresponds 
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FIG. 3: Precessing motion in the vortex front and the vortex 
bundle behind it. (Top) Oscillations in the NMR signal (to- 
tal transverse magnetization M± [3) when the front arrives 
to the pick-up coil result from the precessing motion of vor- 
tices. (Bottom) Fourier transforms of parts of the NMR signal 
marked in the top panel: The front (red) and the cluster be- 
hind it (blue). Since the NMR coil is fixed to the rotating 
frame, the precession frequencies in the laboratory frame are 
counted from the angular velocity f2, as shown by the arrows. 



to the maximum of the power signal, while most of the 
energy is released only after this moment. Thus by the 
time the front has reached the end of the sample, the 
vortex configuration behind the front still includes less 
vortices than in the equilibrium state. The shoulder at 
ilt ^ 1000 and the duration of the signal hint that the 
late relaxation behind the front is the laminar spin-up 
of the superfluid component [l3|. The fit of the experi- 
mental record at 17 — 1.2rad/s to a simple model which 
accounts for the turbulent front creating 0.35 of the equi- 
librium number of vortices, followed by the laminar re- 
laxation with a time constant of 500 s, is shown by the 
green curve in Fig.[2]and is in reasonable agreement with 
the measurement. 

Precessing vortex motion:-Independent confirma- 
tion of the low density of vortices behind the front comes 
from observations of precessing vortex motion around the 
axis of the cylinder, Fig. |31 The precession is visible as 
oscillations of the NMR signal when the vortex arrange- 
ment is not perfectly axially symmetric, which often hap- 
pens in the front and immediately behind it. For vortex 
clusters with uniform density the ratio of their precession 
frequency fig to f2 is equal to the ratio of the number of 
vortices in the cluster to that in the equilibrium state. 
For the front itself the precession frequency is about half 
of fig 13 ■ Thus the measurement in Fig. [3] shows that at 
0.2 Tc the cluster behind the front has about 0.4 of the 
equilibrium vortex density, in agreement with the ther- 
mal signal. 
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FIG. 4: Vortex front in numerical calculations. Snapshots 
of vortex configurations at different temperatures are ana- 
lyzed when the front is located at the same axial position 
z = SAmm. (Top) Azimuthal velocity {va)^ = ilsV as a func- 
tion of z, when averaged over the vortices in the cross sec- 
tion < r < 0.8 R. The vortex-free state is before the 
front {z > 35 mm) and the quasi-equilibrium vortex state 
with r2s(r) < SI is behind the front. (Insert) The val- 
ues of Sis behind the front together with a simple model 
Eq. ^ and the precession frequency of the vortex ends in 
the front region (2 = 30 35 mm) plotted as a function of 
the mutual friction parameter a{T). (Bottom) Turbulence 
in the front is supported by vortex reconnection, which have 
been counted here within bins of width Az — 2 mm over a 
time interval in which the front propagates 0.5 mm. Param- 
eters: cylinder length 40 mm, R — 1.5 mm, S7 — l.Orad/s, 
a(r) = 37.2 exp (-1.97 Tc/T). 



Numerical calculations:-To understand the ob- 
served reduction in the number of vortices numerical sim- 
ulations have been performed using the vortex filament 
model and the Biot-Savart law for computing Vg [3]. 
Above 0.3 Tc, such simulations demonstrate that the av- 
erage vortex number per cross section stays constant as a 
function of z behind the front and corresponds to solid- 
body rotation with the angular velocity of the con- 
Below 0.3 Tc, the vortex density behind the 
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front also stays constant as a function of z, but corre- 
sponds to solid-body rotation with a velocity fig < 
(Fig. 21). In other words, the vortex system develops its 
own rotating frame. The velocity of the vortex frame 
rig decreases prominently with decreasing temperature 
and at 0.20 Tc agrees with the experimentally determined 
value of about 0.4 fi. fJg is found to depend weakly on the 
initial conditions, i.e. on the number and configuration of 
the seed vortices which start the front motion. A similar 
dependence is observed also experimentally. The calcu- 
lations confirm that the precession frequency of the front 
is approximately rJs/2 (insert in Fig. 2]). The difference 
in the rotation velocities of the front and the vortex bun- 
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die behind it is made possible by vortex reconnections, 
as seen in the lower panel of Fig. 01 These reconnections 
lead to the steady-state turbulence in the front. 

Decoupling of rotating frames: The emergence 
of a new rotating frame with angular velocity fig can 
be qualitatively understood on the basis of the coarse- 
grained hydrodynamic equation for Vg, 

^ + V^i - V, X (V X V,) ^ F„ + Fa . (2) 

Unlike classical fluids, two forces appear here on the rhs: 
the mutual friction force between the normal and super- 
fluid components, 

Fa = -aid X ((Vs - Vn) X (V X Vs)) , (3) 

and the force due to the line tension [isj . 

Fa = -A(V X V,) X (V X cj) . (4) 

Here w is a unit vector along the vorticity V x Vg, a is the 
dissipative mutual friction parameter (while the reactive 
mutual friction force with the parameter a' is neglected 
at low T). The line tension parameter A is given by 
A ~ {k / 4tt) hi{£ / a) , where k is the circulation quantum, 
a is the vortex core diameter, and i the inter-vortex dis- 
tance. Two different Reynolds numbers, Rea and Rca, 
are needed to characterize superfluid vortex flow at fixed 
normal component. These parameters represent the rela- 
tive magnitude of the inertial term in Eq. Q with respect 
to the mutual friction (F^) and line tension (Fa) forces: 

Re„ « 1/a , Rca = UR/X = nR^/X . (5) 

The two Reynolds numbers control different regimes of 
the two fluid hydrodynamics: for the onset of turbulent 
flow is required Re a ^ 1 and Re^ ^ 1, while quasi- 
classical Kolmogorov turbulence is found in the regime 
Rcq < (Rca)^^^ and the non-structured quantum tur- 
bulence in the Vincn regime 

The Reynolds numbers describe the two competing 
tendencies experienced by the vortices behind the front: 
spin-up due to mutual friction from the normal compo- 
nent, and spin-down by the precessing vortex front, which 
rotates slower than the vortex bundle behind the front. 
The interplay of these mechanisms leads to the establish- 
ment of a quasi-equilibrium vortex state behind the front. 
This state is in solid-body rotation with the velocity fls, 
which is in general a function of the two Reynolds num- 
bers, changing from the full coupling {fig = fl) at high 
temperatures to the decoupled state (fJs 0) in the 
T — !■ limit. As a simple interpolation between these 
two regimes we may suggest 

a/f^ = (l+Rc„/ReA)~\ (6) 



which is in agreement with the results of the numerical 
simulations (Fig.lH insert). 

Conclusions:-Both measurements and numerical cal- 
culations demonstrate that the vortex state behind 
the propagating vortex front corresponds to quasi- 
equilibrium solid-body rotation at a velocity fig which 
is less than the angular velocity of the normal compo- 
nent: Qs < ^- Thus the angular velocity of the frame in 
which the vortices are at rest deviates from that of the 
container. The decoupling of the two frames becomes 
more and more prominent below 0.3 Tc with increasing 
Reynolds number Re^. This would lead to a total de- 
coupling in the limit T — > 0, unless a different effect in- 
tervenes. Possible candidates for such zero-temperature 
coupling include the circular-motion analog of the Unruh 
effect and the interaction of the vortex-core-bound 

quasiparticle excitations [l7| with the container bound- 
aries. 
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